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We investigate one-dimensional (ID) and 2D radial domain- wall (DW) states in the system of 
two nonlinear-Schrodinger/Gross-Pitaevskii (NLS/GP) equations, which are coupled by the lin- 
ear mixing and by the nonlinear XPM (cross-phase-modulation). The system has straightforward 
applications to two-component Bose-Einstein condensates, and to the bimodal light propagation 
in nonlinear optics. In the former case, the two components represent different hyperfine atomic 
states, while in the latter setting they correspond to orthogonal polarizations of light. Conditions 
guaranteeing the stability of flat continuous wave (CW) asymmetric bimodal states are established, 
followed by the study of families of the corresponding DW patterns. Approximate analytical solu- 
tions for the DWs are found near the point of the symmetry-breaking bifurcation of the CW states. 
An exact DW solution is produced for ratio 3 : 1 of the XPM and SPM coefficients. The DWs be- 
tween flat asymmetric states, which are mirror images to each other, are completely stable, and all 
other species of the DWs, with zero crossings in one or two components, are fully unstable. Interac- 
tions between two DWs are considered too, and an effective potential accounting for the attraction 
between them is derived analytically. Direct simulations demonstrate merger and annihilation of 
the interacting DWs. The analysis is extended for the system including single- and double-peak 
external potentials. Generic solutions for trapped DWs are obtained in a numerical form, and their 
stability is investigated. An exact stable solution is found for the DW trapped by a single-peak 
potential. In the 2D geometry, stable two-component vortices are found, with topological charges 
s = 1,2,3. Radial oscillations of annular DW-shaped pulsons, with s = 0, 1,2, are studied too. A 
linear relation between the period of the oscillations and the mean radius of the DW ring is derived 
analytically. 

PACS numbers: 42.65.Tg; 03.75.Lm; 05.45.Yv; 47.20.Ky 
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I. INTRODUCTION 



A ubiquitous type of topologically protected patterns in binary (two- component) nonlinear systems is represented 
by domain walls (DWs), alias "grain boundaries". A commonly known origin of DWs is in the theory of media with a 
vectorial local order parameter, such as magnetics [l], ferroelectrics 0, and liquid crystals [3 . In systems described by 
binary wave functions (at the fundamental or phenomenological level), the DW represents a transient layer between 
semi-infinite domains carrying different components, or distinct combinations of the components. 

A simple but physically significant example is a rectilinear border between two regions occupied by spontaneously 
emerging roll structures with different orientations. This is, for instance, a generic defect observed in patterns on the 
surface of thermal-convection layers [IHS] . The theoretical description of such patterns is based on systems of coupled 
Ginzburg-Landau equations, each equation governing a slowly varying amplitude of a plane wave whose superposition 
forms the DW [7 . Formally similar systems of coupled nonlinear- Schrodinger (NLS) equations and Gross-Pitaevskii 
(GP) equations describe, respectively, the co-propagation of electromagnetic waves with orthogonal polarizations 
in nonlinear optical fibers [8 , and binary mixtures of Bose-Einstein condensates (BECs) in cigar-shaped traps [9j. 
While the Ginzburg-Landau systems are dissipative ones, on the contrary to the conservative NLS/GP systems, their 
stationary versions essentially coincide, hence DW patterns, generated by the stationary equations, are ubiquitous, 
playing a fundamental role in sundry physical media, both dissipative and conservative ones. Of course, the dynamics 
of perturbed DW structures may be different in the dissipative and conservative systems. 

The basic models which give rise to DWs formed by two fields feature nonlinear XPM (cross-phase-modulation) 
interaction between the fields. In this context, one-dimensional (ID) solutions for DWs were reported in Ref. [7] for 
grain boundaries in thermal-convection patterns, in Refs. [TO^.TT for temporal-domain DWs between electromagnetic 
waves with orthogonal circular polarizations in bimodal optical fibers, and (in a number of different forms) in two- 
component BECs [12 . Various extensions of these settings were investigated too, including grain boundaries between 
domains filled by different cellular [13 and quasi-periodic [14 patterns, and between traveling-wave domains in the 
model of the oscillatory thermal convection [15 . Optical DWs between polarized waves with different wavelengths 
were also studied . Moreover, the analysis was performed for DW states in the discrete version of the nonlinearly- 
coupled NLS system, which describes arrays of parallel bimodal optical waveguides [TT or BEG fragmented in a 
deep optical-lattice potential (in the latter case, the system takes the form of the Bose-Hubbard model) [18]. Also 
considered were the transition to immiscibility in binary BEG with long-range dipole-dipole interactions [19 , and 
various forms of DWs in the three-component spinor BEG [2Q| . 

As concerns the experiment, linear grain boundaries between patches filled by rolls with different orientations had 
been well documented in many observations of the thermal convection [5 , 21 . Similar linear defects were reported in 
laser cavities [22 and in experimental studies of BEG [S^. Well-pronounced DW structures have also been created 
in bimodal optical fibers [23j and in fiber lasers [24]. 

In many physically relevant settings, the nonlinear interaction between coexisting waves competes with the linear 
interconversion between them, which strongly affects the DW patterns in such binary systems. In particular, the linear 
coupling between the orthogonal polarizations induced by the twist or elliptic deformation of the optical fiber gives 
rise to an effective force accelerating the corresponding DW [11] . The linear interconversion between two components 
of a binary BEG representing different atomic hyperfine states, coupled by a resonant radio- frequency wave, affects 
the formation of the DW in immiscible binary condensates [26]. The objective of the present work is to develop a 
comprehensive study of DW states in systems of NLS/GP equations coupled by both the linear and nonlinear (XPM) 
terms, in ID and 2D geometries. Related two-component vortical states in 2D are studied too. 



The paper is organized as follows. The model is formulated in Section jllj In Section |III| we consider flat continuous- 
wave (GW) states, which represent backgrounds supporting DWs. Both symmetric and asymmetric GW states are 
found, and the symmetry-breaking bifurcation (SBB), which gives rise to asymmetric GW backgrounds, is identified. 
The stability of the flat states is also investigated in this section. Basic types of DWs in the free space (without 



an external potential) are considered in Section IV, where both analytical and numerical solutions for the DWs are 
reported, and the stability of the DW patterns is established. Section |V] deals with one and two DWs interacting with 
an external potential, in the form of one or two peaks. In particular, an analytical stable solution is produced for a 



DW pinned to a potential peak. Two-dimensional axisymmetric patterns are considered in Section VI Stable vortices 



with topological charges s = 1,2,3, supported by the asymmetric GW background, are found, and radial oscillations 



of annular pulsons in the form of circular DWs are studied too. The paper is concluded by Section VII 
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II. THE MODEL 

Our starting point is the system of scaled one- dimensional NLS/GP equations for two wave functions ?/^i^2, coupled 
by the linear and nonlinear (XPM) terms: 



i = - (1/2) {i,^)^^ + a IVi 1^1+5 IV'2|' V'l - «V'2, 



i{^2\ = -(1/2) (^2),, + <T|t/'2|'V'2+5|V'l|'V'2-«V'l, 



(1) 



where n is the rate of the linear interconversion between the two atomic states, if the system is interpreted in terms 
of BEG [26 . The same equations, with time t replaced by the propagation distance, and x replaced by the reduced 
time, r = t — z/Vgr, where Vgr is the group velocity of the carrier wave, may be realized in optics as the model of 
the light propagation in an ordinary or photonic-crystal fiber [8]. In the fiber-optic model, amplitudes ?/^i^2 represent 
two mutually orthogonal polarizations of light, with the linear interconversion induced by the birefringence or twist 
of the fiber, for the circular or linear polarizations, respectively. Coefficients a and g in Eqs. ([T]) account for the SPM 
(self-phase- modulation) and XPM nonlinearities, respectively (in the ordinary optical fiber, the XPM/SPM ratio is 
g/a = 2 and 4/3 for the of circular and linear polarizations, respectively). 

To secure the modulational stability of CW states supporting DWs, coefficient g will be kept positive, which 
corresponds to the repulsive XPM nonlinearity, while SPM coefficient a may have any sign. In optics, opposite signs 
of the XPM and SPM coefficients is an exotic situation, which is, nevertheless, possible in photonic-crystal fibers [8j. 
In BEG, the sign of either coefficient may be switched by means of the Feshbach-resonance effect [9 . 

The value of g may be fixed by dint of an obvious rescaling [for instance, it is possible to set ^ = 3, which is a 



natural choice in view of the existence of an exact DW solution in the form of Eqs. (36) and (37), see below, which 
requires g = 3a]. In addition, the rescaling allows one to fix = 1, thus we will assume hz = ±1. Nevertheless, in 
the analytical expressions written below, we keep g and K: as free parameters, as it is easier to analyze the results in 
such a form. 

Stationary solutions to Eqs. ([T]) with chemical potential /i are sought for as 

i^i,2{x,t)=e-"''<P,,2{x), (2) 
with real functions 4'i,2{x) satisfying equations 

IJ.4>i + (1/2) 4>'l - a(Pl - g4>l(Pi + K(P2 = 0, 

(3) 

li<j)2 + (1/2) <j)'^ - <j4>1 - g(t>l4>2 + k4>i = 0, 

with the prime standing for d/dx. The energy (Hamiltonian) corresponding to stationary states (jij) is iJ = T-Ldx, 
with density 

H = (i/2) [((/>;)' + ((/>^2f 

+ (a/2) + 4) + 9^l4>l - '2i^M2- (4) 



III. FLAT (CONTINUOUS- WAVE) STATES 

Symmetric flat (x-independent) solutions to Eqs. ([3|, with equal amplitudes of both components, are 

(l)i=(j)2 = Ao = ^{ii^K) / {(j^g). (5) 
The Hamiltonian density (|4| for this solution is 

'^^symm = - / (^ + Cr) . (6) 

Antisymmetric CW solutions are equivalent to Eq. ([5| with k replaced by —k. In view of this relation, we will define 
the symmetric solutions as those for n = +1, while the antisymmetric ones will be replaced by symmetric states for 
n = -1. 
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Asymmetric fiat states, which are generated by the SBB, can also be found in the exact form: 



(7) 



2(7 y4(T2 (5 -a)' 
with the signs of (j>i and (j>2 determined by relation 

<Pi4>2 = K{g- cr)"-^ , (8) 
which is compatible with Eqs. The Hamiltonian density ^ for these states is 

n^symm = mV (2(t) - / (q - a) . (9) 

Taking into account the above convention, <pi(p2 > 0, and the condition that <p1 2 must be positive, it follows from 
Eqs. ^ that the asymmetric CW state emerges, i.e., the SBB takes place, at a specific value of the chemical potential, 

H = 2aK{g - = licT, (10) 
the corresponding value of Aq at the bifurcation point being obtained by the substitution of this value into Eq. 

Al = K{g-a)-\ (11) 
The asymmetric solution ([8|, ^ exists at /i^ > /i^^. Slightly above the bifurcation point, i.e., at 

M = Mcr + (^M, with < |/icr| , (12) 

one can expand Eqs. ([7| and (|8| in powers of small 5/i, which yields 



^1,2 



" ±^a/^ + 0(<5m). (13) 



g — cr 2 V cr 



A. The modulational stability of the symmetric solution at the bifurcation point 

A crucial condition necessary for the existence of stable DWs is the absence of the modulational instability (MI) of 
the corresponding CW background. Here we explicitly consider the MI of the symmetric CW, and, in particular, we 
will find conditions providing for the stability of the symmetric state exactly at the SBB point. Examining this case 
secures that the asymmetric CW states are not subject to the MI - at least, close enough to the bifurcation point. 

To analyze the MI, we look for perturbed solutions to Eqs. ([T]) in the well-known general form, i.e., as 

i^i,2{x,t) = [Ao ^ ai^2{x,t)]exp[-i/j.t ^ ixi,2 {x,t))] , (14) 

where amplitude Aq is the same as in Eq. (I5|, while ai^2 and Xi,2 are infinitesimal perturbations of the amplitudes and 
phases of the two components (this form of the solution implies the incorporation of the small complex perturbations 
into the unperturbed CW state). Then, eigenmodes of the perturbations are sought for as 

{ai,2(x,t),xi,2(^,^)} = |a^°^, Xi°2}exp (7^ + ^P^), (15) 

where p is an arbitrary real wavenumber of the perturbation, and 7 the corresponding instability gain, which may be 
complex. The stability condition is that Ke{j{p)} = for all real p. 

The substitution of expressions ([l4| and ([T5| into Eqs. ([T]) and linearization with respect to infinitesimal pertur- 
bations yields the following system of equations, which actually splits into two separate subsystems for ^a^^^ + <^2^^^ 7 



[2^2 ^ ^) ^ (i/2)p^] (a^ + ) + 7^0 (x^ + ) = 0; 



(16) 
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0, 



[2A^ ig-a)-2n- (1/2) p^] (a^ - ) - 7^0 (x^ " ) = 0. 



The resolvability conditions for Eqs. (16) and (iTTfc yield, respectively, the following expressions for j{p): 



7i = -(l/2)p2 [2{g + a)Al + {l/2)p'], 
7! = (2/t + (1/2) p2) [2 {g -a)Al-2K- (1/2) p^] . 
The stability condition ensuing from Eq. ( [l8| , i.e., 7^ < 0, is obvious: 

g + a> 0. 

Expression (19) simplifies at the bifurcation point (10), where Aq = k/ {g — a), as per Eq. (11): 



7- 



■(l/2)p2 [2K+(l/2)p2 



(17) 

(18) 
(19) 

(20) 

(21) 



Evidently, the stability condition following from Eq. ([21| is > 0, which, as a matter of fact, means that only the 
symmetric flat solution may be stable at the SBB point, while its antisymmetric counterpart, that (as defined above) 
corresponds to < 0, is unstable. Further, from > and the positiveness of expression (11) for A^^, condition 
g — a > follows. Combined with Eq. (20), it gives rise to the following relation between the XPM and SPM 
coefficients necessary for the existence and stability of the asymmetric CW states: 



(22) 



while <j may be positive or negative. In fact, it will be demonstrated below that DW solutions do not exists for a < 0. 
Note that condition (22) does not hold for ^ = 0, hence the DWs that can be found in the system with the solely 



linear coupling are always unstable, as demonstrated in Ref. |27j. In that work, it was shown that the unstable DW 
gives rise to an expanding layer filled with turbulent waves. 



B. The modulational stability of the asymmetric background 



In the general case of the asymmetric CW background, the perturbed solution is taken as [cf. Eq. (14)] 

i^iM^^ t) = Ai^2 [1 + ^i,2(^, t)] exp [-iiit + zxi,2 {x, t))] , (23) 
where Ai^2 are given by Eqs. ([7|), and infinitesimal perturbations are taken as [cf. Eqs. (15)] 

{6i,2(x, t), xi,2(^, t)} = {6^% Xi'2} exp (7^ + ipx) . (24) 



The substitution of expressions ( [23| ) and (24) into Eqs. ([!]), and the subsequent linearization, lead to a system of four 
linear equations: 



7^2 



(0) 

7X1 



(0) 

7X2 



(AO) 

'A, 




= 0, 






= 0, 






= 0, 


^1 


-5f) 


= 0. 



(25) 



In the special case of Ai = A2 = Aq, it is easy to check that Eqs. (25) are tantamount to Eqs. (16) and (17) written 
above. 

The dispersion relation between 7 and is determined by the resolvability condition of system (25): 



7 








7 







7 




= 0. 



(26) 
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Equation ( [26| ) was solved numerically, yielding four roots 7(p^). It has been verified that, with Ai and A2 taken as 
per Eqs. ^ and ([8|, and for all > 0, all the roots satisfy condition Re {7 (p^) } =0, if k, is positive and inequality 
(22) holds. Thus, as long as the symmetric CW are stable at the SBB point, its asymmetric counterparts are stable 
too. 



IV. DOMAIN- WALL SOLUTIONS IN THE FREE SPACE 



Numerical results 



Stationary solutions to equations ([3| were constructed by means of the Newton-Raphson method for the corre- 
sponding nonlinear boundary-value problem. In particular, the boundary conditions were fixed as zero values of the 
derivatives at both edges of the integration domain. 

DWs are built as transient layers fusing together CW states of different types. Obviously, such patterns are possible 



at /i^ > /i^^ [see Eq. (10)], where asymmetric CW solutions exist, as given by Eqs. ([7|. Figure [l| shows that, in this 



case, one can find four different types of the transient layers, if dark solitons are counted too. These types may be 
classified by values of coupled fields (^1, ^2) in the uniform states connected by the DWs: 



{(01 (X = -00), 02 {X -00)) , (01 (X = +00), 02 {X 

{(Ai, A2) , (A2, Ai)} ; {(Ai, A2) , (-A2, -A,)} ; 
{(Ai, A2) , (-Ai, -A2)} ; {(Ao, Ao) , (-Ao, -Ao)} 



-00))} 



(27) 



types in Eq. (27) [Figs. 



[recall that Ap is given b y Eq. ([5|, and Ai^2 are given by Eqs. ([t])] . In f act, only the patterns of the first and second 

1(a) and |(b)| are true DWs, while Figs. |l(c)| and |(d)[ corresponding to the profiles of the 
third and fourth types in Eq. ([27|), are paired dark solitons. It is also clear that the dark-soliton pairs of the latter 



type, displayed in Fig. l(d)[ are unstable past the SBB point, as symmetric CW state ([5| is unstable in this case. 

It is relevant to mention that DWs connecting the symmetric and asymmetric states - for instance, (^1,^2) and 
(Aq^Aq) - are impossible, because a stationary DW may only exist between two asymptotic fiat states with equal 
Hamiltonian densities [11 . Comparing the respective densities (|6| and ([9|, one can immediately conclude that they 
coincide solely at the bifurcation point (10). 



To investigate the stability of the DW patterns, small perturbations were added to the stationary solutions: 

0l(x, t) = 0i(x) + vi{x)e-'^' + <(x)e^^*^ 
0^(x, t) = 02(x) + V2{x)e-'^' + ^5(x)e^^*^ 



(28) 



where ^i, Ui and ^2, U2 constitute eigenmodes of the infinitesimal perturbation, and 7 is the corresponding eigen- 
frequency that, in general, may be complex. Substituting expressions (28) into Eqs. (|3| and linearizing around the 
stationary solutions leads to the following eigenvalue problem. 



Li 

^0201 
-^0201 4 



^0102 - ^ 
-^0102 
-L2 



^0102 
#102 - 
(70^ 
L9 













Ui 




Ui 




h 










\U2 / 




\U2) 



(29) 



£2 = n+(l/2)(f/dx^ 



9<t>l 



(30) 



This eigenvalue problem can be solved using a simple finite-difference scheme. Accordingly, the solution is identified 
as a stable one if all the eigenfrequencies are real. 

The stability analysis outlined above demonstrates that the DW family of the first type in Eq. (27), which is 
represented by Fig. |l(a)[ is completely stable (past the SBB point, where it exists, along with the asymmetric CW 
states). In this case, all the other types of the DW and dark-soliton solutions [i.e., all of them which cross zero at least 
in one component, see Fig. [TJb,c,d)] are completely unstable^ due to the presence of imaginary eigenfrequencies in the 
spectrum of small perturbations. We stress that, unlike the trivial background instability of the pattern of the last 
type in Eq. (27), those of the second and third types, which are represented by Figs. |l(b)| and |(c)[ are destabilized by 
perturbations localized around the transient layer, while the CW background is stable. Direct simulations demonstrate 
that these unstable DWs decay into expanding turbulent patterns (see Fig. |2|. The stability of the DW of the first 
type in Eq. (27) was also verified by direct simulations (not shown here). 
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FIG. 1. Typical examples of the domain walls (a, b) and paired dark solitons (c, d), found for^ = 3, (7 = ^^=1, and /x = 2. 
The four panels represent the patterns of the four types defined in Eq. (27). Stable and unstable solutions are depicted by the 
continues and dashed lines, respectively. 



B. An approximate solution for the DW near the bifurcation point 



An analytical solution for the DW can be constructed in an approximate asymptotic form near the bifurcation 
point (10), (11), i.e., for ja taken as in Eq. (12). In this case, the approximate solution can be sought for as 

^l(x) = Acr + S(l)i{x) + S(l)2{x), 

M^) = Act - SM^) + H2{X), (31) 

where it is implied that S(j)i ~ ^/SJl and S(j)2 ^ Sfi, cf. Eq. (13). Substituting expressions (31) into Eqs. ^ and 
expanding the result in powers of Sjii yields a relation between S<pi and 5(j)2 at order Sja, 



2{h + k) 

Next, at order the expansion yields the equation for S(pi{x): 



(32) 



4- Sfi • S(l)i 



16a {d^if = 0. 



(33) 



An exact solution of the DW type to Eq. (33) is 



6(l)i{x) 



1 Sfi 



tanh 



2- — —djix ) . 



(34) 
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FIG. 2. (Color online) Examples for the evolution of u nstab le patterns of the second (a) and third (b) types from Eq. (27) 
The initial conditions and the parameters are as in Fig. 1(b) and|(c)| for (a) and (b), respectively. 



With regard to the CW-stability condition (22), we conclude that solution (34) exists if the sign of Sji coincides 



with the sign of a, and only for ^/x > 0, which means that the solution exists solely for a > 0. In fact, the numerical 
solution demonstrates too that DWs cannot be found at a < 0, even if asymmetric CW states ([t]) exist in this case, 
with /i < 0. 

In what follows below, we will also need an expression for the full density of the DW solution. The above formulas 
yield 



[<P^{x)f + [<P2{x)f 
2k, Sji g — (J Sfi 
(J q -\- (J (J 



sech 



a 



dfix 



(35) 



Examples of the analytically predicted profiles, as given by Eqs. (31), (32) and (34), together with their numerically 
found counterparts, are displayed in Fig. |3(a)^ for parameters ^ = 3, <j = 1,/>: = 1 [i.e., /icr = 1, see Eq. (10)] and 
^/x = 0.1, 0.5, 1 and 2. For the same examples, the difference between the numerical and approximate results is shown 
in Fig. |3(b)[ The results presented in Fig. 3(a) demonstrate that the prediction loses its accuracy with the increase 
of Sjii. On the other hand, comparing the core of the analytical and numerical solutions (the transient layer), one can 
see that the analytical results are not necessarily most accurate for small Sjii. Actually, this approximation is most 
suitable for intermediate values of J/i. 



C. The exact solution for g — ?>(t 

A particular exact solution to Eqs. (|3| for the DW can be found by means of the following ansatz: 

(j)i{x) = Uo -\-Ui tanh (Ax) , 

(36) 

02 (^) = Uq — Ui tanh (Xx) . 

The substitution of the ansatz into Eqs. ([3| demonstrates that it yields an exact solution at g = 3<j, for chemical 
potential /i = + , with coefficients 

A^ = /i - = (4cr)"^ {2k + A^) = (4cr)"^ (/i + , 

= (4a)"' A' = (4a)"' (/i - k) (37) 
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FIG. 3. (a) Comparison between the analytical approximation given by Eq. (34), which is shown by dotted lines, and 
numerically found profiles of the domain walls (solid lines), for g = 3, a=l, k^I (/Xcr = 1) and S/j = 0.1, 0.5, 1 and 2. (b) 
The differences between the numerical and approximate analytical results, A(/)i,2 = (0i,2)numer — (0i,2)approx, for the examples 
presented in panel (a). 



(under condition /i — > 0). In fact, this particular exact stationary solution is similar to the one that was found, also 
in the exact form, for the stationary version of coupled Ginzburg-Landau equations (but without the linear coupling) 
in Ref. [7 . Note that relation g = 3cr, which is necessary for the existence of the exact solution, complies with the 
CW-stability condition (22). In terms of this exact solution, /i may be considered as a free parameter, i.e., we actually 
have a one-parameter solution family, under constraint g = 3a. An example of the exact DW solution can be seen 
in Fig. lfa)[ where the parameters were chosen to comply with the exact solution. Finally, the total local density of 
solution (36) is 



<Pl{x) 



2 [{U^ + Ul) - U^sech^ (Xx)] , 



(38) 



cf. Eq. (35). 



It may also be relevant to mention that, under condition g = 3<7, Eqs. (|3| has another formal DW-like solution, 
given by Eqs. ( 36 ) and ( 37 ) with tanh replaced by coth. However, this additional solution is singular, therefore it has 



no physical meaning. 



V. DOMAIN WALLS IN THE PRESENCE OF EXTERNAL POTENTIALS 

The model based on Eqs. ^ can be naturally extended to include an external potential, W{x): 



ficPi + (1/2) 0'/ - - g<Pl<Pi + K<P2 - W{x) 



+ (1/2) (p2 - ^'Pl - 9<Pl<p2 + K^i - W{x) 



0, 



(39) 



In this section we report the analysis performed for the DW patterns supported by this setting, with W{x) representing 
single- or double-peak potentials, with the intention to predict pinning of the DWs by such potentials. Stable pinning 
by potential maxima (rather than minima) may be possible because the DWs core features a minimum of the total 
density, see Eqs. (35) and (38), hence the total energy of the system may be minimized by placing the core around a 
local maximum of the potential. 
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A. The exact solution in the presence of a single-peak potential 



It is possible to find an exact solution to Eqs. (|39| for the DW if the potential is chosen as 

W{x) = Wo sech^ (Ax) 



(40) 



where Wq and A are considered as given parameters. The corresponding solution can be looked for in the form of the 
same ansatz (36) as used for finding the exact solution in the free-space setting. The substitution of the ansatz into 



Eqs. (39) demonstrates that it yields an exact solution in the present case under the following condition imposed on 



parameters of the system: 



3A^ + 2Wq 
A2 + 21^0 



(41) 



which goes over into the above relation, g = 3a, in the limit of Wq = 0. Further, coefficients Uq, Ui, and /i of the 
exact solution are given by the following expressions: 



A-2 {\^ + 2Wo) (A^ + k) 
(4<tA2)-' (A2 + 2I^o) (a" 
(4(t)-^ {\^ + 2Wo). 



2k] 



(42) 



It is easy to see that, in the limit of Wq = 0, Eqs. (|42| carry over into the above exact solution, given by Eqs. ( |37[ ). 
Thus, this exact solution is an extension of the previous one, although it has no free parameters [note that A, which 
was an adjustable parameter of the free-space solution, is now fixed by the given form of potential ([4Q|)]. 

As argued above, the present exact solution is expected to be stable if W{x) represents a repulsive potential barrier 
(peak), with Wq > 0, and the solution should be unstable, against spontaneous escape from the pinned state, in the 
case of the attractive potential well, with Wq < [cf. Eq. ( [48| ) below]. In turn, Eq. (41) with T^o > gives g < 3a. 
Note also that, unlike its counterpart in the free space, this exact solution may exist at a < 0: In the case of the 
potential weh {Wq < 0), with A^ < —2Wo < 3A^, Eq. (41) yields a < 0. However, this solution should be unstable 
according to the above argument. 



B. Analysis of the interaction between DWs and the interaction of the DW with the external potential 

The interaction of the DW with an external potential can be investigated in an approximate form. In fact, a similar 
problem which also admits an approximate analytical treatment is the interaction between two broadly separated 
DWs with opposite polarities (i.e., mirror images of each other) in the free space, therefore we start with this case. 

Assuming that the two DWs are set at distance L which is large in comparison with the inner width of each DW, 
the interaction between them can be analyzed by means of the method elaborated in Ref. To this end, the 

approximate expression for the nearly fiat fields in the region between the far separated solitons is taken as 

^n{x) = A, + 2Un [exp (-2A \x - 6|) + exp (-2A \x - 6|)] , (43) 

where n = 1, 2, and ^41^2 are given by Eqs. (I7|), ^1,2 are coordinates of the centers of the two DWs, so that L = ^2 — ^i, 
and the exponential terms represent small decaying tails of the DWs on top of the fiat background. The decay rate 
A > can be found in the general case, but the expression for it is cumbersome; amplitude Ui^2 are not known in an 
exact form in the general case, as they may only be found from full solutions for individual DWs. In the special case 



of the exact solution given by Eqs. (36), at ^ = 3cr), coefficients A and Ui^2 in Eq. (43) are actually given by Eq. (37): 



A = v7^^, Ui = U2 = {l/2)^{fi-i^) /a. (44) 

Then, using Hamiltonian produced by density (|4| and the method developed in Ref. [28 (identifying the term ac- 
counting for the interaction energy in the expression for the full Hamiltonian) , the effective potential of the interaction 
between the two separated DWs is found in the following form: 

Uint{L) = -8A {Ui + /7|) exp (-2AL) . (45) 
In the case of the exact DW solution given by Eqs. ([36| and (|37|, this expression takes an explicit form, 

Uint {L) = -4 (/i - exp {-2^/jr^L) . (46) 
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A simple but essential property of expressions (45) and (46) is that they obviously predict attraction between the two 
DWs. 

With the external potential W{x) included into Eqs. (39), Hamiltonian density (j4| is modified as 



n 



(0'l) +(02 



+grKPl - 2^^102 + 2W{x) + <pl-A\- Al) , 



(47) 



where constants A\ 2 are subtracted from 0^ 2 (^) for convenience (to cancel a formally diverging constant term in the 
Hamiltonian), ^41^2 being the same as in Eqs. This means that the energy of the interaction of the DW with the 
external potential is 



pot 



+00 



W{x) [<Pl{x) + <Pl{x)-A\-Al]dx. 



(48) 



Assuming that W{x) represents a broad potential barrier or well with a width much larger than the thickness of 
the DW, and treating the external potential as a perturbation (i.e., neglecting the distortion of the solution under the 
action of the potential), the substitution of the exact DW solution given by Eqs. (36) and (37) into Eq. (48) readily 
yields: 



(49) 



where ^ is the coordinate of the center of the DW. A straightforward consequence of Eq. (49) is that, as conjectured 
above, the DW tends to be trapped at local maxima of the external potential W{x), as, due to sign minus in Eq. 
( 49 ) , they correspond to minima of the effective potential ( 49 ) . 

If two DWs are trapped at two particular maxima of W(a;"i and ^2, separated by large distance L, the equilibrium 
condition for each DW is the vanishing of the total force produced by the interaction of the DW with its counterpart 
and with the external potential: 



d_ 
W2 



2^- 



x3/2 



W (C2) + 4 (Ai - nf^ exp (-2V^7^ (C2 - 6)) 



exp (C2 - 6)) 



0. 



(50) 



For example, if the potential is periodic, W{x) = ecos (27rx/A), with large period A, one may consider the pair of 
DWs trapped at two adjacent potential maxima. Substituting this potential into Eqs. (50), it is easy to check that its 
minimum strength, necessary for holding the DW pair (i.e., preventing it from merger due to the mutual attraction) 
is Cmin = (2/7r) crA (/i - f^f^'^ cxp {-^/JjT^K) . 



C. Numerical results for the model with an external potential 



The above analysis predicted that two DWs created in the free space attract each other. Direct simulations confirm 
the prediction, see a typical example in Fig.[4j Eventually, the two DWs annihilate into the stable uniform asymmetric 
state. Naturally, for larger initial values of the separation between the DWs, the attraction force is weaker, and a 
considerably longer time is required for the DW pair to manifest the interaction. 

In the presence of the external potential in Eqs. (39), the numerical investigation was carried out for two different 
shapes of W{x). First, 



the single potential barrier was considered, taken in the same form ( |4Q| ) which was used to 
obtain the exact solution. The corresponding stationary DW solutions are similar to their counterparts found in the 
free space. Specifically, out of the four types of the DWs categorized by the values of ((/)i(x), ^2(^)) at x = ±00 as per 
Eq. (27), only the first one is stable, see a typical example in Fig. [sja). In addition to the DWs, the potential barrier 
supports states in the form of "bubbles" (defined as per Ref. [29]), i.e., patterns wit h identica l values of (0i(x), 02(^)) 
at X = ±00, and localized perturbations of ^i,2(^) around the barrier, see Figs. 5(b) -(c) Naturally, the bubbles 



supported by the asymmetric and symmetric CWs are stable and unstable, respectively. 

The expectation that the DW is stably trapped by the potential peak is confirmed by direct simulations displayed 
in Fig. |6(a)[ for the DW of the first type, in terms of Eq. (27), i.e. {(^41742), {A2A1)} [it is the same DW which is 
displayed in Fig. [sja)]. The DW, if shifted from the potential maximum by Ax = 5, performs decaying oscillations 
around the peak. The decay of the oscillations is explained by emission of radiation waves into the background by 
the oscillating DW. 
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FIG. 4. The simulated evolution of a pair of identical domain walls, initially separated by distance 2Ax — 10, is displayed by 
means of density contour plots, for g = 3, a — 1, n — 1, and ii — 2. 
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FIG. 5. Typical profiles of the domain-wall (a) and "bubble" (b,c) stationary solutions found in the presence of the potential 
peak ( |4Q| , for g — 8 /3, a = 1 , — 1, A = 1, Wo = 0.1, and /x = 2.4 [these values of the parameters admit the existence of the 
exact DW solution (41)-(42), which is actually displayed in (a)]. The stable and unstable solutions are depicted by continuous 
and dashed curves, respectively. 



The attraction of DWs to the potential peak is further illustrated by Fig. (sjb), which displays a pair of identical 
DWs symmetrically placed, at t = 0, at distances Ax = ±5 from the peak. Due to the attraction of both DWs to the 
peak, in this case their annihilation happens much sooner than for the same pair in the free space, cf. Fig. |4]for the 
same value of separation 2 Ax between the DWs. 

The numerical analysis was also carried out for the double-peak potential, taken in the form of 

W{x) = Wo {2x/Lf exp (l - {2x/Lf^ , (51) 
where is the height of the two peaks and L is the distance between them. The numerical solution of Eqs. (I39|) 



with potential (51) has revealed multiple structures, which may be considered as transient layers between CW states 
filling the three regions separated by the two potential barriers. These structures can be categorized into two groups, 
built of symmetric or asymmetric CWs. The structures of the former type are arranged as the following sets of the 
CW states in the three regions [cf. Eq. ([27[)]: 

{(Ao, Ao) , (Ao, Ao) , (Ao, Ao)} ; {(Aq, Aq) , (Aq, Aq) , (-Aq, -Aq)} ; 

{(Ao, Ao) , (-Ao, -Ao) , (Ao, Ao)} . (52) 



These three varieties may be interpreted, respectively, as containing none, one, or two dark solitons trapped in each 
component. 
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FIG. 6. (a) The evolution of the domain wah shifted by Ax = 5 from the maximum of potential barrier (40), for the same 
parameters as in Fig.jsja). (b) Similar to panel (a), but for the pair of DWs, placed symmetrically on both sides of the potential 
maximum, at x — ±5. 



Using various combinations of the asymmetric CW states, we have built the following patterns supported by the 
double-peak potential, cf. Eqs. ([27|) and (52): 



{{A,, A2) , {A,,A2) , (Ai, A2)} ; {(^1,^2) , (Ai, A2) , {A2, A,)} ; 
{{A^,A2),{A2,A^),{A,,A2)}; 

{(^1,^2) , (Ai, A2) , {-AU-A2)} ; {(^1,^2) , {-AU-A2) , (Ai, A2)} ; 
{(^1, ^2) , {Ai,A2) , {-A2, -A,)} ; {(-^1, -A2) , iAuA2) , {A2, Ai)} ; 
{{Ai,A2) , {-AU-A2) , (^2, ^1)} ; {{Ai,A2) , iA2, ^1) , i-Au -A2)} ; 
{{A,,A2),{-A2,-A^),{AuA2)}. 



(53) 



(54) 



The three patterns (53) feature, respectively, none, one, or two trapped DWs, while seven patterns (54) include dark 



solitons or zero-crossing (sign-changing) DWs. As before, those patterns which do not cross zero in any component 
are stable [in the case of the first pattern in (52), this is, naturally, true prior to the SBB], while all the solutions 
featuring at least one zero crossing are unstable. Examples of the patterns of these types are presented in Fig. [7| for 
Wo = 0.1 and L = 10. 

The attractive interaction of the DWs with the potential peaks is additionally illustrated by simulations displayed 
in Fig. |8] In particular, if a single DW is initially placed at the midpoint between the peaks, the evolution of this 
obvious ly un stable configuration leads to damped oscillations of the DW around either peak to w hich i t is pulled, 
see F ig. |8(a)[ On the other hand, symmetrically placing pairs of identical DWs either between (Fig. |8(b)[ ) or outside 
(Fig. |8(c)[ ) of the potential peaks, we observe coherent oscillations with the apparently unbroken symmetry. In the 
latter case, the symmetric oscillations feature very slow damping, which is explained by the small value of l^o = 0.1 in 
this case. Similar simulations for larger Wq demonstrate oscillations which converge much faster to a stable symmetric 
configuration, with each DW trapped by one potential peak. 
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FIG. 7. Typical examples of patterns supported by the double-barrier potential ([51|, with Wo = 0.1 and L = 10, constants ^, a, 
K, and /i being the same as in Fig. [4] Panels (a), (b) and (c) demonstrate the stable patterns of the following types, in terms of 
Eq. ([53}: {(Ai, As) , (A2, Ai) , (Ai, A2)}, {(Ai, A2) , (Ai, A2) , (A2, Ai)}, and {(Ai, A2) , (Ai, A2) , (Ai, A2)}, respectively. The 
unstable pattern of type {(^o, ^o) , (^o, ^o) , (^o, ^0)}, from set ([52|, is presented in panel (d). 



VI. VORTICES AND TWO-DIMENSIONAL CIRCULAR DOMAIN WALLS 



The 2D version of model ([T]) is 

i{^i)t = - {1/2) [(V^i), 

i{^2)t = - (1/2) [{^2)a 



^V^2, 



(55) 



In terms of BEG, Eqs. ([55| admit the straightforward interpretation as the GP equations for the two-component 
condensate in a 2D pancake-shaped configuration. In terms of optics, these equations, with t replaced by propagation 
coordinate govern the transmission of a stationary beam through the bulk nonlinear medium, with functions 
'01,2 representing two circular polarizations. In the latter case, the linear mixing between the polarizations can be 
induced by the linear birefringence, which, in turn, may be imposed by mechanical stress applied perpendicular to 
the propagation axis |8J, or, alternatively, by dc magnetic field applied in the same direction |30J (the birefringence 
imposed by the transverse magnetic field leads to the classical Gotton-Mouton/Voigt effects [3T]). 

General axisymmetric solutions to Eqs. (55) are looked for in the form of the optical vortices (alias 2D dark solitons) 
[HI [32] , which also correspond to the vortex modes in BEG [9 : 



V^i,2(x,^,t) = 0i,2(r,t)exp(i56>)exp(-z/it), 



(56) 



where r and are the polar coordinates in the (x, y) plane, and integer s is the topological charge (vorticity, alias 
"spin"). 



15 






(a) 



(b) 



(c) 



FIG. 8. The evolution of a single DW and pairs of DWs in the presence of the double-peak external potential, (51), with 
Wo = 0.1 and L = 10. The other parameters are as in Figs. |4] and |7| In panel (a), the single DW is initially positioned exactly 
between the two peaks. The evolution of pairs of DWs, symmetrically placed at x = ±2.5 and x = ±7.5, is displayed in panels 
(b) and (c), respectively. 



The substitution of expressions (56) into Eqs. (55) leads to the equations for radial wave functions 0i,2 



i{4>i)t + ^l4>l + (1/2) [{4>i)rr + r-\4>i)r - sV-20i] 

-(j{4>lf -g{4>2?(t)l+K(t)2 = 0, 
i{4>2)t + IJ-4>2 + (1/2) [{4>2)rr + r"^(02)r " S^'^^a] 

-(j{4>2? -g{4>if(t)2 + >^(t>i = 0. 



(57) 



Obviously, stationary CW solutions in the 2D model are identical to those obtained for the ID model, see Eqs. ([5| 
and In particular, the bifurcation described by Eqs. (10) and (11) is relevant in the 2D case too. 



The study of the existence, stability and dynamics of 2D axisymmetric patterns, generated by Eqs. (57), is presented 
below in two steps. First, stationary vortices, supported by the asymmetric CW background, are obtained for spins 
s = 0, 1, 2 and 3, and their stability is determined. Then, the evolution of the pulsons^ i.e., circular DWs oscillating in 
the radial direction, is investigated by means of direct simulations. In this connection, it is relevant to mention that 
patterns in the form of circular DWs are well known in various magnetic media [33] 



A. Stationary vortices 



To find stationary solutions of Eqs. (57), we proceed to the time-independent version of these equations and apply 
the Newton-Raphson method to the respective nonlinear boundary-value problem. The boundary conditions demand 
that 01,2(0 ^ at r ^ 0, while, at r ^ oo, 0i,2(^) asymptotically approach the symmetric or asymmetric CW 
solutions, ([5| or ([t]). Examples of the radial profiles of such modes are displayed in Fig.joJ For the s = 0, no stationary 
solutions exist, apart from the obvious fiat states (which are not shown in Fig. |9|. For each non-zero value of the 
spin that was examined, 5 = 1,2 and 3, two families of nontrivial stationary solutions were obtained, supported [past 
bifurcation point (10)] by the symmetric and asymmetric CW states at r ^ oo, hence they demonstrate precisely the 
same existence and bifurcation features as the corresponding fiat states, ([5| and ([7|. 
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FIG. 9. (a,b,c): Radial profiles of the symmetric and asymmetric stationary vortices, with spin s = 1, 2, and 3, respectively 
(no non-flat stationary 2D solutions were found for s = 0). The parameters are g = 3, a = 1, k = 1, and ii — 2. Solid and 
dashed curves depict, respectively, two components of the stable asymmetric solutions, and unstable symmetric ones. 

To explore the stability of the 2D stationary solutions, we seek for perturbed solutions as [cf. Eqs. 



(28])] 



^i(r,6>,t) = U{r) + [/7+(r)exp(m6>) + /7_(r) exp(-m6>)] exp(7nt), 
Mr,0,t) = V{r) + [F+(r) exp(ml9) + F_(r) exp(-ml9)] exp(7^t), 



(58) 



where integer n > is an arbitrary azimuthal index of the perturbation, and 7„ is the corresponding instabihty growth 
rate. After substituting expressions (58) into Eqs. (jSTl) and hnearizing, the following system of equations is obtained: 



fiU+ + ij„U+ + {1/2)U'; + (2r)-' U: 

-gUV{Vl + V+) - aU'^iUl + 2U+ 
HU- + i-fnU- + {l/2)U'l + (2r)-^ U: 

-gUV{V* + V-)- aU^iUl + 2U- 
fiV+ + i7„y+ + (1/2)1^;' + (2r)-' V: 

-gUV{Ul + U+) - aV^{Vl + 2V+ 
HV- + ij„V- + {1/2)V!! + (2r)"' V: 

-gUV{Ul + U-) - crV^{V* + 2V. 



-{s + nf{2r^)-'U+ 

■ gV'^U+ + kV+ = 0; 
-(s-n)2(2r2)"V_ 

■ gV^U- + kV- = 0; 
-{s + nf (2r2)-V+ 

■ gU^U+ + kU+ = 0; 
-(s-n)2(2r2)-V_ 
- gU'^V- + kU- = 0, 



(59) 



where the prime stands for d/dr. We treat Eqs. (59) as an algebraic eigenvalue problem for 7^ and solve it directly, 
using a finite-difference method. The largest instability- growth rate is identified as the real part of the most unstable 
eigenvalue, max{Re(7n)}. Following this approach, we have confirmed that the only source of the destabilization is 
the transition between the symmetric and asymmetric modes, which, as said above, is actually driven by the SBB 
in the fiat background at r ^ 00. In particular, no azimuthal instability, that would break the axial symmetry, 
was found for s = 1,2,3, for all integer values of n considered. The absence of the azimuthal instability in the case 
of the self-defocusing nonlinearity is actually a natural feature [34]. In addition, this finding implies that the fiat 
(quasi-one-dimensional) DW is stable against corrugations in the 2D setting. 



B. Direct simulations — shrinking domain walls and pulsons 



Another aspect of the 2D model based on equations ( 55 ) was examined by investigating the development of circular 



DWs into pulsons, i.e., annular grain boundaries periodically shrinking and expanding in the radial directions. It is 
known that, while exact solutions for such pulsons do not exist, in some models, such as the 2D single-component 
sine-Gordon equation, the pulsons may be remarkably stable, featuring hundreds [35]|37j and thousands [36 of radial 
pulsations with very little loss. 
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FIG. 10. The initial evolution of the pulsons (circular domain walls) for s = (a), s = 1 (b), 
of the circular DW is i?o = 20, the other parameters being g = 3, a=l, k=1 and /i = 2. 



and 



2 fc). The initial radius 



We have performed simulations of pulsating radial DWs for three values of the spin, 5 = 0, 1 and 2. The simula- 
tions were implemented by means of the linearized Crank-Nicolson scheme. As initial conditions, the following DW 
configuration, based on the asymmetric uniform states ([7|), was used. 



h{r,t 



0) = (1/2) tanh"(r) {Ai [1 + tanh(r - i^o)] + ^2[1 - tanh(r - Rq)]} . 
0) = (1/2) tanh^(r) {Ai[l - tanh(r - i?o)] + ^2[1 + tanh(r - Rq)]} . 



(60) 



where Rq is the initial radius of the DW. Note that this initial ansatz complies with the necessary boundary condition 
at r ^ 0, 01,2 (^) ~ The boundary conditions at the right edge of the integration interval, < r < p, were 
adopted as {d(j)i^2/d^ + d4'i,2/dt) \r=p = 0. These conditions prevent the reflection of the emitted radiation from the 
boundary. 

For each value of the spin considered here, s = 0,1,2, the simulations were run for several values of the initial 
radius, Rq. Typical examples of the initial stage of the observed evolution are demonstrated in Fig. [To) for Rq = 20 
and for the parameters g = 3, a = l^hz = l and /i = 2. Similar results were obtained for other values of Rq. 

It was observed that, while the pulsons shrink and expand in a quasi-periodic manner, their smallest and largest 
radii, i^min and i^max, do not remain constant, slowly decreasing from a cycle to a cycle, as shown in Fig. 11 In 
particular, the decay of i^max is roughly exponential in time, with some irregularities observed at t 400 and at the 
final stage of the evolution. The decay is plausibly caused by the emission of radiation waves by the pulsating DW. 

Also decreasing is the oscillation period, r^, taken as the time interval between two consecutive points at which the 
pulson expands to the largest radius, ri = t{r = i^max,/) — ^(^ = ^max,^-!), where / is the number of the cycle. The 
relation between the slowly decreasing period and the effective maximum radius, i^max,eff = (^max,/ + ^max,^-i)/2, 



is presented in Fig. 12 The plots demonstrate a nearly linear dependence for all s. The same linear dependence was 
obtained for different values of the initial radius, including Rq = 10, 15,30, and 40. 

The linear relation between ri and i^max can be easily explained. Indeed, considering a large-radius circular DW, 
with radius R much larger than the thickness of the DW in the radial direction, one can define the effective mass for 
the radial pulsations, M = 27: Rm^ and the effective surface-tension energy, £^st = 27: Ra^ where m and a are the 
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FIG. 11. The temporal decay of the largest and smallest radii of the pulson, in the cases shown in Fig. |1Q| for s = 0, 1 and 2. 
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FIG. 12. The relation between the gradually decreasing oscillation period of the pulson, and the effective maximum radius, 
i?max,eff, for s = 0, 1 and 2. The parameters are the same as in Fig. IIOI 



effective mass and surface-energy densities of the quasi-one-dimensional DW. Thus, the Newton's equation of motion 
for the DW in the radial direction is 



d / dR 
— 27rmR-- 
dt V dt 



_d^ 
dR 



{27TaR) , 



(61) 



from where the law of motion follows: R'^{t) = i^^ax — {(^/m)t'^^ assuming that the motion starts with the zero initial 
velocity and R = i^max- According to this result, the shrinking DW ring will bounce back from the center at the 
moment of time ri/2 = ^/a/mRmax^ which is obviously equal to a half of the period. This result explains the linear 
proportionality of ri to i^max- 

Finally, Fig. [13] displays the final stage of the evolution, observed when the pulson's radius and oscillation period 
have been sufficiently reduced. In this case, the oscillatory behavior fades o ut and the pulsons transform into the 
stable asymmetric stationary states which are presented above in subsection VI A| (in particular, they are flat for 
s = 0, and feature the vortex shape for s = 1 and 2). For the case considered above (with Rq = 20), the transition to 
the eventual stationary state occurs at t ~ 830, for all the values of the spin. 
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VII. CONCLUSIONS 



The objective of this work was to investigate fundamental DW ( domain- wah) modes in the system of nonhnear- 
Schrodinger/Gross-Pitaevskii equations, coupled by the linear and XPM (cross-phase-modulation) terms. The system 
of coupled equations has natural realizations in two-component Bose-Einstein condensate, and in nonlinear optics, 
where the two wave functions represent amplitudes of co-propagating signals with orthogonal polarizations. First, 
conditions providing for the stability of the uniform CW (continuous wave) symmetric and asymmetric bimodal states, 
which support the DW patterns, were identified, and then the general families of DW solutions were constructed in 
the ID setting. In particular, approximate analytical solutions were found near the symmetry-breaking bifurcation 
point of the CW states, and an exact solution was found for the XPM/SPM ratio 3:1. The DW states connecting 
asymptotically fiat asymmetric states which are mirror images to each other (without the change of the overall sign) 
are stable, while all other types of the DWs, which feature zero crossings (including dark solitons), are unstable. 
Interactions between two DWs with opposite polarities were also considered. The potential of the attraction between 
them was found in the analytical form, and numerical simulations have demonstrated that the attraction leads to 
annihilation of the DW pair. 

DWs trapped by the single or double potential peaks were investigated too. An exact solution for the DW trapped 
by a single peak was found. In the general case, it was predicted and corroborated by direct simulations that the 
bound state of the DW placed at a maximum of the external potential is stable. The interaction of two DWs in the 
presence of the single- and double-peak potentials was also studied. 

The analysis was extended to axisymmetric patterns in the 2D geometry, including vortices carrying the topological 
charge s = 1,2,3 and supported by the asymmetric fiat background at infinity. Stable stationary states for the 
vorticity-carrying DW rings were found. The stability of the vortices against azimuthal perturbations was verified 
through the computation of the corresponding eigenvalues. Oscillations of annular DW-shaped pulsons were studied 
by means of systematic simulations (the linear relation between the period of the radial oscillations and largest radius 
of the annular DW was obtained in an analytical form). The evolution of the pulsons ends up with their relaxation 
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into the stationary 2D modes (vortices, for 5 7^ 0, or simply the fiat asymmetric state in the case of 5 = 0). 

This work may be extended by a more systematic investigation of the 2D system, without assuming the axial 
symmetry of the patterns. In particular, oscillations of eccentricity in elliptically deformed pulsons may be interesting 
to study, cf. Ref. [37 . In ID, it may be also interesting to study in detail the behavior of DWs against the backdrop 
of periodic potentials (optical lattices), as well as in nonstationary systems with time-dependent parameters. 



ACKNOWLEDGMENT 

The work of J.Z. was supported, in a part, by a postdoctoral fellowship provided by the Tel Aviv University, and 
by grant No. 149/2006 from the German-Israel Foundation. 



[1] V. G. Bar'yakhtar, M. V. Chetkin, B. A. Ivanov, and S. N. Gadetskii, Dynamics of Topological Magnetic Solitons (Springer- 

Verlag: Berlin, 1994). 
[2] D. Damjanovic, Rep. Progr. Phys. 61, 1267 (1998). 

[3] P. G. de Gennes and J. Prost, The Physics of Liquid Crystals (Oxford University Press: New York, 1995). 

[4] P. Manneville and Y. Pomeau, Phil. Mag. A 48, 607 (1983); G. Tesauro and M. C. Cross, ibid. 56, 703 (1987); T. Passot 

and A. C. Newell, Physica D 74, 301 (1994). 
[5] E. Bodenschatz, W. Pesch, and G. Ahlers, Ann. Rev. Fluid Mech. 32, 709 (2000). 

[6] A. A. Nepomnyashchy, M. G. Velarde, and P. Colinet, Interfacial Phenomena and Convection (Chapman and Hall/CRC: 
New York, 2001). 

[7] B. A. Malomed, A. A. Nepomnyashchy, and M. I. Tribelsky, Phys. Rev. A 42, 7244 (1990). 

[8] Y. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic Press: San Diego, 2003). 

[9] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Clarendon Press: Oxford, 2003). 
[10] M. Haelterman and A. P. Sheppard, Phys. Lett. A 185, 265 (1994). 
[11] B. A. Malomed, Phys. Rev. E 50, 1565 (1994). 

[12] M. Trippenbach, K. Coral, K. Rzazewski, B. Malomed, and Y. B. Band, J. Phys. B: At. Mol. Opt. 33, 4017 (2000); P. 
Ohberg and L. Santos, Phys. Rev. Lett. 86, 2918 (2001); S. Coen and M. Haelterman, ibid. 87, 140401 (2001); J. J. 
Garcia-Ripoll, V. M. Perez-Carcfa and F. Sols, Phys. Rev. A 66, 021602 (2002); P. G. Kevrekidis, H. E. Nistazakis, D. 
J. Frantzeskakis, B. A. Malomed, and R. Carretero- Gonzalez, Eur. Phys. J. D 28, 181 (2004); B. A. Malomed, H. E. 
Nistazakis, D. J. Frantzeskakis, and P. G. Kevrekidis, Phys. Rev. A 70, 043616 (2004); J. Ruostekoski, ibid. 70, 041601 (R) 
(2004); K. Kasamatsu and M. Tsubota, Phys. Rev. Lett. 93, 100402 (2004); K. Kasamatsu, M. Tsubota, and M. Ueda, 
ibid. 93, 250406 (2004); P. G. Kevrekidis, H. Susanto, R. Carretero-Gonzalez, B. A. Malomed, and D. J. Frantzeskakis, 
Phys. Rev. E 72, 066604 (2005); Y.-C. Kuo, W.-W. Lin, and S.-F. Shieh, Physica D 211, 211 (2005). 

[13] C. Kubstrup, H. Herrero, and C. Perez-Carcia, Phys. Rev. E 54, 1560 (1996); Y. N. Young and H. Riecke, Physica D 163, 
166 (2002). 

[14] H. Rotstein and B. A. Malomed, Physica Scripta 62, 164 (2000). 

[15] B. A. Malomed, Physica Scripta 57, 115 (1997); H. Sakaguchi and B. A. Malomed, Physica D 118, 250 (1998). 
[16] S. Wabnitz, IEEE Phot. Tech. Lett. 21, 875 (2009). 

[17] P. G. Kevrekidis, B. A. Malomed, D. J. Frantzeskakis, and A. R. Bishop, Phys. Rev. E 67, 036614 (2003); Z. Rapti, A. 

Trombettoni, P. G. Kevrekidis, D. J. Frantzeskakis, and B. A. Malomed, and A. R. Bishop, Phys. Lett. A 330, 95 (2004). 
[18] A. Mering and M. Fleischhauer, Phys. Rev. A 81, 011603(R) (2010). 

[19] G. Gligoric, A. Maluckov, M. Stepic, L. Hadzievski, and B. A. Malomed, Phys. Rev. A 82, 033624 (2010). 

[20] W. Zhang, D. L. Zhou, M.-S. Chang, M. S. Chapman, and L. You, Phys. Rev. Lett. 95, 180403 (2005); M. Uchiyama, 

J. leda, and M. Wadati, J. Phys. Soc. Jpn. 75, 064002 (2006); H. Saito, Y. Kawaguchi, and M. Ueda, Phys. Rev. A 75, 

013621 (2007); H. E. Nistazakis, D. J. Frantzeskakis, P. G. Kevrekidis, B. A. Malomed, R. Carretero-Gonzalez, and A. R. 

Bishop, ibid. A 76, 063603 (2007); Z.-D. Li, Q.-Y. Li, P.-B. He, J.-Q. Liang, W. M. Liu, and G. Fu, Phys. Rev. A 81, 

015602 (2010). 

[21] V. Steinberg, G. Ahlers, and D. S. Cannell, Physica Scripta 32, 534 (1985); Y. C. Hu, R. Ecke, and G. Ahlers, Phys. Rev. 
E 48, 4399 (1993). 

[22] Ye. Larionova, U. Peschel, A. Esteban-Martin, J. Garcia Monreal, and C. O. Weiss, Phys. Rev. A 69, 033803 (2004). 
[23] S. Pitois, G. Millot, and S. Wabnitz, Phys. Rev. Lett. 81, 1409 (1998). 

[24] H. Zhang, D. Y. Tang, L. M. Zhao, and X. Wu, Phys. Rev. B 80, 052302 (2009); H. Zhang, D. Y. Tang, and L. M. Zhao, 
Opt. Exp. 18, 4428 (2010). 

[25] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattor, and D. M. Stamper-Kurn, Nature 443, 312 (2006); N. Gemelke, 

X. Zhang, C.-L. Hung, and C. Chin, ibid. 460, 995 (2009). 
[26] D. T. Son and M. A. Stephanov, Phys. Rev. A 65, 063621 (2002); B. Deconinck, P. G. Kevrekidis, H. E. Nistazakis, and 

D. J. Frantzeskakis, ibid. A 70, 063605 (2004); M. L Merhasin, B. A. Malomed, and R. Driben, J. Phys. B: At. Mol. Opt. 

38, 877 (2005); A. Niederberger, B. A. Malomed, and M. Lewenstein, Phys. Rev. A 82, 043622 (2010). 



21 



[27] R. S. Tasgal and B. A. Malomed, Physica Scripta 60, 418 (1999). 

[28] B. A. Malomed and A. A. Nepomnyashchy, Europhys. Lett. 27, 649 (1994); B. A. Malomed, Phys. Rev. E 58, 7928 (1998). 
[29] I. V. Barashenkov and E. Y. Panova, Physica D 69, 114 (1993); I. V. Barashenkov, Phys. Rev. Lett. 77, 1193 (1996). 
[30] J. F. Dillon, J. P. Remeika, and C. R. Staton, J. AppL Phys. 41, 4613 (1970); Y. Linzon, K. A. Rutkowska, B. A. Malomed, 

and R. Morandotti, Phys. Rev. Lett. 103, 053902 (2009); K. A. Rutkowska, B. A. Malomed, and R. Morandotti, Opt. 

Exp. 18, 8879 (2010). 
[31] R. W. Boyd, Nonlinear Optics (Academic Press, San Diego: 2008). 

[32] P. Coullet, L. Gil, and F. Rocca, Opt. Commun. 73, 403 (1989); G. A. Swartzlander and C. T. Law, Phys. Rev. Lett. 69, 
2503 (1992); L V. Basistiy, V. Y. Bazhenov, M. S. Soskin, and M. V. Vasnetsov, Opt. Commun. 103, 422 (1993); D. Rozas, 
C. T. Law, and G. A. Swartzlander, J. Opt. Soc. Am. B 14, 3054 (1997); M. S. Soskin, V. N. Gorshkov, M. V. Vasnetsov, 
J. T. Malos, and N. R. Heckenberg, Phys. Rev. A 56, 4064 (1997); D. M. Palacios, L D. Maleev, A. S. Marathay, and G. 
A. Swartzlander, Phys. Rev. Lett. 92, 143905 (2004). 

[33] S. A. Langer, R. E. Goldstein, and D. P. Jackson, Phys. Rev. A 46, 4894 (1992); K. O. Ng and D. Vanderbilt, Phys. Rev. 
B 52, 2177 (1995); R. P. Cowburn, D. K. Koltsov, A. O. Adeyeye, M. E. Welland, and D. M. Tricker, Phys. Rev. Lett. 83, 
1042 (1999). 

[34] B. A. Malomed, D. Mihalache, F. Wise, and L. Torner,. J. Optics B: Quant. Semicl. Opt. 7, R53 (2005). 

[35] L L. Bogolyubskii and V. G. Makhankov, JETP Lett. 25, 107 (1977); P. L. Christiansen, P. S. Lomdahl, and N. J. Zabusky, 

Appl. Phys. Lett. 39, 170 (1981); E. M. Maslov, Phys. Lett. A 151, 47 (1990); A. G. Bratsos, J. Comp. Appl. Math. 206, 

251 (2007). 

[36] J. Geicke, Physica Scripta 29, 431 (1984). 

[37] P. L. Christiansen, N. Gr0nbech- Jensen, P. S. Lomdahl, and B. A. Malomed, Physica Scripta 55, 131 (1997). 



